Far-infrared spectroscopy of spin excitations and Dzyaloshinskii-Moriya interactions 
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We have studied spin excitation spectra in the Shastry-Sutherland model compound SrCu2 (603)2 
in magnetic fields using far-infrared Fourier spectroscopy. The transitions from the ground singlet 
state to the triplet state at 24cm~^ and to several bound triplet states are induced by the electric 
field component of the far- infrared light. To explain the light absorption in the spin system we 
invoke a dynamic Dzyaloshinskii-Moriya (DM) mechanism where light couples to a phonon mode, 
allowing the DM interaction. Two optical phonons couple light to the singlet to triplet transition 
in SrCu2(B03)2. One is a-polarized and creates an intra-dimer dynamic DM along the c axis. The 
other is c-polarized and creates an intra-dimer dynamic DM interaction, it is in the {ab) plane and 
perpendicular to the dimer axis. Singlet levels at 21.5 and 28.6 cm~^ anti-cross with the first triplet 
as is seen in far-infrared spectra. We used a cluster of two dimers with a periodic boundary condition 
to perform a model calculation with scaled intra- and inter-dimer exchange interactions. Two static 
DM interactions are sufficient to describe the observed triplet state spectra. The static inter-dimer 
DM in the c-direction di = 0.7 cm~^ splits the triplet state sub-levels in zero field [Cepas et al., 
Phys. Rev. Lett. 87, 167205 (2001)]. The static intra-dimer DM in the (ab) plane (perpendicular 
to the dimer axis) ci2 = 1.8 cm~^ , allowed by the buckling of CuBOs planes, couples the triplet state 
to the 28.6 cm~^ singlet as is seen from the avoided crossing. 

PACS numbers: 75.10.Jm, 78.30.Hv, 71.70.Gm, 76.30.Fc 



I. INTRODUCTION 

In spin systems with a ground singlet state and ex- 
cited triplet state the energy gap between the singlet 
and the triplet can be tuned with an external magnetic 
field. In SrCu2 (603)2 it was discovered that in mag- 
netic fields above 22T, where the spin gap is expected to 
close, several magnetization plateaus appear,^ At magne- 
tization plateaus the triplets form a pattern which breaks 
the translational symmetry of the crystal structurerSi The 
heavy mass of the triplet excitations arising from an al- 
most flat dispersion of energy on momentum'^ favors the 
build-up of magnetic superstructures. Below the critical 
field SrCu2 (603)2 has a ground state described first by 
Shastry and Sutherland]— 

SrCu2 (603)2 consists of planes of CU6O3 and Sr 
atoms between the planes. Cu^+ spins (S=l/2) form 
Cu-Cu dimers arranged into orthogonal dimer net- 
work. SrCu2 (603)2 is an experimental realization of 
a Shastry-Sutherland modeli^ In the model there is an 
anti-ferromagnetic intra-dimer exchange coupling ji and 
inter-dimer coupling j2 between spins on the nearest- 
neighbor dimers (Fig.^. In the limit of a = J2/J1 = 
the problem reduces to that of isolated dimers where 
the ground state is the product of singlet states and the 
first excited triplet state is at energy At = ji above the 
ground state, where At is the energy per dimer. Shas- 



try and Sutherland showed that for < a < 0.5 sin- 
glets on all dimers is an exact ground state too. The 
exactness of the ground state and the heavy mass of 
triplet excitations is the consequence of frustration orig- 
inating from the special geometry of the dimer lattice 
in the Shastry-Sutherland model where the bonds on 
neighboring dimers are orthogonal. Later on it has been 
shown that singlets on all dimers is the exact ground 
state for a larger range of a up to the quantum critical 
point ac ~ 0.7. At the quantum critical point the spin 
gap vanishes and a long-range anti-ferromagnetic order is 
established. Different theoretical approaches have been 
used to calculate ofc (see Refi^ for review) . It is possible 
that between the exact singlet ground state and the anti- 
ferromagnetic state in certain range of a other gapped 
spin states exist— iSiSiiSiii 

The singlet-triplet gap in SrCu2 (603)2, At — 
24cm~^, has been measured directly by several exper- 
imental techniques: inelastic neutron scatteringSiiS, elec- 
tron spin resonanceiSii^ (ESR), Raman scatteringi^, and 
far-infrared (FIR) spectroscopy:^^ Additional informa- 
tion besides At is needed to determine the exchange pa- 
rameters of SrCu2 (603)2. The dispersion of the triplet 
excitation is not informative because of its flatness^., but 
positions of other excited states or the temperature de- 
pendence of thermodynamic parameters can be used for 
determining the exchange parameters. Miyahara and 
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UedaA found ji = 59cm~^ and a = 0.635. They added 
an interlayer coupling = 0.09ji to the model to obtain 
a better fit of the magnetization T- dependence above the 
critical temperature ksT > At- Based on the analysis 
of excitation spectra^^ ji — 50cm~^ and a — 0.603 were 
proposed. Such scattering of parameters could be either 
due to the incomplete model or due to the approxima- 
tions made in theoretical calculations. SrCu2 (603)2 is 
near to the quantum critical point ac where the energy 
levels of the spin system are sensitive to the choice of ji 
and a. A singlet level in the spin gap at 21cm~^ found 
in the ESR spectraM may help to find proper parameters 
for the model. 

Interactions other than inter- and intra-dimer ex- 
change coupling can spoil the exactness of the ground 
state. This is important in high magnetic fields where 
the triplet state becomes degenerate with the ground 
singlet state. At this critical field even a weak inter- 
action between the singlet and the triplet state mixes 
the two states completely. The singlet and triplet 
state anti-crossing effects were seen in the high field 
ESR experiments.^^ A possible anti-symmetric interac- 
tion which couples the singlet and the triplet states is 
the DM interaction. An intra-dimer DM is allowed by 
symmetry but its strength is not known below room 
temperature. Above room temperature d2 = 2.5 cm~^ 
has been estimated from the ESR linewidthii^ The inter- 
dimer DM interaction, di = 1.5 cm~^, perpendicular to 
the dimer planes^^ partially lifts the degeneracy of the 
triplet state but does not couple the triplet state to the 
singlet state. The effect of DM interactions on the mag- 
netic dipole active ESR transitions in SrCu2 (603)2 was 
investigated theoretically in Ref.^^. 

Lattice distortions, static or dynamic, are important 
in SrCu2 (603)2 since they lower the crystal symmetry 
and allow magnetic interactions which are otherwise for- 
bidden in a more symmetric environment. SrCu2 (603)2 
has a structural phase transition at 395 K^- that induces 
a buckling of CU6O3 planes in the low T phase. As 
the phase transition point is approached from below the 
Raman-active 62cm~^ optical phonon mode softensiSi 
Acoustic phonon modes have spin-phonon coupling at 
magnetization plateaus;^^: ft has been proposed that a 
spin superstructure at 1/8 plateau observed by nuclear 
magnetic resonance at 35 mK is stabilized by a lattice 
distortion.^ Instantaneous breaking of lattice symmetry 
by an optical phonon allows electric dipole active singlet- 
triplet transitionsS^ that explains FIR polarized absorp- 
tion spectra in SrCu2 (603)2 fi^- 

Our aim is to find out which additional interactions 
are required to the Shastry-Sutherland model that add 
triplet corrections to the ground state. For that we 
do FIR absorption measurements with polarized light 
in magnetic field and compare the absorption line fre- 
quencies and intensities with values calculated with a two 
dimer model including the dynamic DM effect. The im- 
portant information is in the polarization and magnetic 
field dependence of the FIR absorption lines and in the 



TABLE L Singlet and triplet excitations observed in the FIR 
spectra at 4.4 K in the order of increasing zero field energies 
hujo (in cm~^ units). When a line is visible in two Ei polariza- 
tions, both are indicated. The corresponding Hi polarizations 
are also indicated. S and T label the singlet and triplet states; 
+ ( — ) denotes levels which energy increases (decreases) with 
Bo and indicates levels where the energy stays constant; Qa 
and Qc are the 3- factors with Bo || a and Bo |[ c, respec- 
tively. The labeling of To levels is shown in Fig.|S] The zero 
field intensities Aq (in cm~^ units) of To (★) are described in 
the text and in figures |1| |21 and |S| High energy excitations 
are labelled by their energies. 
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±0.16 
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0.2 ± 0.1 


42.7 
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±0.16 






T43.5(±) 
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44.7 


±0.4 
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c 
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0.04 ± 0.02 
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2.27 
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a 


b,c 


86 ± 14 


52.24 
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^53. 5 
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b,c 


24 ± 3 


53.44 


±0.07 







avoided crossing effects. 

We studied single crystals of SrCu2(603)2, Ref.— . 
The first sample consisted of two pieces 0.65 mm thick 
in a-direction with the total area of 12 mm^ in the (ac) 
plane. The second sample was 0.6 mm thick in the c- 
direction and had an area of 11.5 mm^ in the {ab) plane. 
The experimental details are described in ReL^^'^^. 



II. RESULTS AND DISCUSSION 
A. FIR spectra and electric dipole transitions 

As the result of the polarization sensitive measure- 
ment of FIR spectra we have identified that the main 
resonances in the spectra are electric dipole transitions, 
rather than being magnetic dipole transitions. In Fig. [21 
differential absorption spectra at 4.4 K relative to 15 K, 
are displayed. The strong absorption lines at 52.3 and 
53.5 cm~^ were identifiedi& as electric dipole transitions, 
that are active in Ei || a polarization. We see the same 
for the 43.0 cm~^ singlet and Tq and Ti triplets (see 
TableU) at 24.2 and 37.5 cm"^, respectively, which are 
present in the spectra measured with Ei || a regardless 
of Hi being perpendicular to the c axis or parallel to it. 
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The lines are missing in Ei || c polarization^-, instead 
a new line appears at 25.5 cm^^, which is identified as 
another component of the triplet Tq. 

The triplets are split by the magnetic field Bq. Differ- 
ential absorption spectra in Ei || a polarization for one 
magnetic field direction, Bq || a, measured relative to the 
zero field, are displayed in Fig. |3| We see an anti-crossing 
of the rom(— ) level with the singlet at 21.5 cm~^ and 
an anti-crossing of the Top{+) level with the singlet S2 
at 28.6 cm~^. All the peaks in the measured spectra in 
different light polarizations and Bp directions were fitted 
with Lorentzians. The results are summarized in Tabled 
and displayed in Figures ^ [SJ and E| The states above 
38cm~^ are labelled by their zero field frequencies. The 
magnetic field independent energy levels are labelled as 
singlets with the exception of those in the middle of the 
triplet levels r(±). 

B. Dynamic Dzyaloshinskii-Moriya mechanism and 
optical transitions: two dimers 

The hamiltonian for a spin pair with exchange coupling 
j and DM interaction d on the bond connecting spins k 
and I reads: 

H^L = (j - 4^)Sfc • + ISfe . dd . Si 

+ d-[SkxSi]+gfiBBo-{Sk + Si). (1) 

Here we included Shekhtman corrections^^'^^ which are 
quadratic in d (see alsc^^). The last term is the Zeeman 
energy of spins in the magnetic field Bq where g is the 
electron spin g-factor and /is is the Bohr magneton. 

The formalism to introduce the spin-phonon coupling 
is similar to one used in Refi^iSi. We are interested in 
singlet to triplet transitions. Therefore the relevant term 
is the anti-symmetric DM interaction d{Q) ■ [Sk x S;] 
which couples the singlet to the triplet state. We expand 
the DM vector d(Q) into a power series of the lattice 
normal coordinate Q 

d(g)=d(0) + |^|Q^oQ + ...., (2) 

where d(0) = d is the static DM interaction in Q. We 
keep terms linear in Q. The full hamiltonian for a spin 
pair including the phonons is 

H'^' = H^i,t + hujp a) a + q{a) + a) dg • [S,. x S;], (3) 

where dg = Iq=o- '^^'^ lattice normal coordinate Q 
is presented in terms of phonon creation and annihilation 
operators and a, Q = q{a'^ +a), where q is the transfor- 
mation coefficient and ujp is the phonon frequency. The 
spin-phonon coupling term in ^ is linear in and a. 
Therefore the phonon states with the occupation num- 
bers n and n' are coupled where n' = n ± 1. We will 
consider only two phonon states |0) and which is jus- 
tified when fc^T <C htUp. 



The normal coordinate Q in the dynamic DM singlet to 
triplet optical transition mechanism belongs to an optical 
phonon. Electric dipole coupling between a phonon and 
light in the long wavelength limit is 

V = eQEi = eq{a^ + a)Ei, (4) 

where e is the effective charge associated with the lat- 
tice normal coordinate Q. Here we assumed Ei || Q and 
dropped the time dependence of V. Once the eigenstates 
of (PI are known the optical transition probability be- 
tween the ground state \(f>) and the excited state |(/)') is 
calculated as /=|(^!)'|V|(/))|2. 

To calculate optical transitions in SrCu2(B03)2 we use 
a two dimer model depicted in Fig.^ In this model intra- 
dimer and inter-dimer superexchange interactions ji and 
j2 are considered. The inter-dimer static DM vector di 
is along the c axis and alternates from bond to bond. 
The intra-dimer static DM vector d2 exists due to the 
buckling of Cu-O-B planesi^ The direction of DM vec- 
tors is defined by the right hand rule where the path is 
along the Cu-O-Cu bond (for di Cu-O-B-O-Cu) in the 
direction of increasing spin index k. In the vector prod- 
uct Sfc X S; the spin with a smaller index is on the left, 
k < I. When a periodic boundary condition is applied 
to the two dimer cluster, bounded by a box drawn with 
a thin dashed line in Fig.^a), an effective spin model is 
obtained where the inter-dimer interactions are doubled, 
Fig.^b). The doubling is necessary to conserve the num- 
ber of next-nearest-neighbor bonds, which is four. 

The hamiltonian for the two dimer cluster is the sum 
of pairwise interactions Q where the sum runs over all 
the bonds in the cluster. We will use a basis \ ABn) where 
A runs over the singlet S and three triplet components 
T^,To, and T-|_ on the ji bond of the dimer (1,2) and 
B over the singlet and triplet states of the dimer (3,4). 
n is the number of phonons, or 1. The basis has 32 
components. Below we consider a- and c-axis phonons, 
shown in Fig.[7| named by the direction of their electric 
dipole moment. 

1. Energy levels 

The effect of the dynamic DM interaction on the po- 
sition of energy levels is small because we take Uojp = 
100 cm~^ that is substantially larger than the singlet- 
triplet gap. We use this value since there are no optical 
phonons with substantial spectral weight below 100 cm^^ 
as our transmission measurements show. The energy 
spectrum can be analyzed separately from the dynamic 
DM effect because of the high phonon energy. The calcu- 
lated energy levels are the same in Fig.2|and|Sl In these 
figures only the zero phonon levels of the triplet Tq and 
S2 are shown. The levels with one excited phonon are 
off-set by fkup to higher energies and are not shown. 

In a two dimer system two singlets, two triplets, and 
a quintet are present. The ground state is a product of 
singlets \SS). The first triplet is a linear combination 
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of \ST) and ITS'). In the two dimer model the singlet- 
triplet splitting is not renormalized by the inter-dimer 
coupling j2 and the energy of the triplet excitation is 
Eto ~ ji- The second singlet, a bound state of two 
triplets, is at Esi = 2ji — 2(2j2)- To stress the fact 
that in the two dimer model with a periodic boundary 
condition the inter-dimer bonds are effectively doubled, 
we write 2j2 explicitly. There are two other bound states 
of two triplets, a triplet at Eti = 2ji — {2j2) and a 
quintet at Eq = 2ji + (2j2). These energies and the 
ground state wavefunction are slightly changed by the 
static DM interactions di and d2- The spin states \STi) 
and \TiS) are strongly mixed by the inter-dimer di since 
they are degenerate in any field. 

The states are labelled in Fig.El The following param- 
eters were used to fit the energy spectra plotted in Fig.^ 
and 13 The energy of one-triplet sublevels rom(O) and 
rop(O) gives us ji = 24.0 cm~^. To get the singlet S2 
at 28.6 cm~^ we use 2j2 = 9.8 cm~^. Triplet levels are 
split in zero field by 2di = 1.4 cm^^. The intra-dimer 
^2 — 1.8 cm^^ induces an avoided crossing of Top(-I-) and 
^2. In a simplified picture the one-triplet excitation is 
the \ST) (or \TS) ) state and the excited singlet is \TT). 
d2 "flips" the singlet to the triplet state on one of the 
dimers and thus couples Top{+) to 82- 



2. c-axis phonon 

The optical c-axis phonon bends the Cu-O-Cu bond 
in the c-direction. We assume that the bending action 
of the phonon is the same on both dimers, Fig.d As 
a result the dynamic DM interaction on the dimer (1,2) 
is Qcdq^ = — (— c?3c,0,0) and on the dimer (3,4) 
dsc = (0, d^c, 0); the orientation of the Cartesian coordi- 
nates is the same as in Fig.n(b). The calculated and the 
measured transition probabilities as a function of mag- 
netic field are plotted in Fig.01(b, d) for two field orien- 
tations. In zero field a line at 25.5 cm~^ is present. The 
area of this line is the only scaling parameter between the 
theory and the experiment. Note that the transition to 
the triplet level, which anti-crosses with 52, is optically 
active when Bq || c. When Bq || a there is no crossing 
for the optically active triplet level. 

The overall agreement between the theory and the ex- 
periment is good. There is a disagreement between the 
intensities of the middle and lower triplet components 
in the theory and in the experiment, Fig.^d). In the 
theory the intensity of the middle component is approx- 
imately three times as strong as the lower component 
while in the experiment they are equal. We tried several 
changes in our model to make the intensities of the two 
triplet components more equal and none of them helped. 
These unfruitful changes were the shift of the phonon 
frequency, a small out-of-plane component of Bq and an 
in-plane component of the inter-dimer DM vector di. 



3. a-axis phonon 

The optical a-axis phonon bends the Cu-O-Cu bond in 
the a-direction and creates a dynamic DM interaction in 
the c-direction, Fig.[3 If we choose Ei || a the dynamic 
DM interaction is created on dimer (1,2), qa^Qa = dsa — 
(0, 0, d^a)- In general, for an arbitrary orientation of Ei 
in the {ah) plane, both dimers will acquire a certain ds^. 
For the time being we assume Ei || a. 

In zero magnetic field the transition to the central 
triplet component is observed, Fig.|21 As the Bq || c field 
is turned on, Fig.|3b, the central line, being a sum of 
two overlapping transitions, conserves its intensity. The 
experimentally observed drop in intensity with increas- 
ing field is a T effect. At 1.8 K (18 T field) the inten- 
sity is recovered. Besides the strong central line there 
are in zero field two sidepeaks ten times weaker at 22.7 
and 25.5 cm~^ corresponding to transitions to the twice 
degenerate states Tom(±) and Top(±). The dynamic 
DM interactions due to the a axis and c-axis phonons 
in this Bq orientation give zero intensity for the side- 
peaks. The detailed analysis of the mechanism causing 
these weak transitions is difficult because in other polar- 
izations and field orientations stronger mechanisms are 
prevailing. The sidepeaks split in the magnetic field and 
an avoided crossing with Si and S2 is seen in the exper- 
iment. 

When the magnetic field is in the (ah) plane two cases 
must be considered, Bq || Ei and Bq -L Ei. In Fig.[31(c, 
d) the Bo II El case is shown. Here are optically active 
the triplet levels which anti-cross with the singlet states. 
In Bq -L El field orientation, Fig.|Hl the optically active 
triplet levels do not anti-cross with the singlet states. 
The mutual orientation of Bq and Ei is important be- 
cause El II a creates daa on the dimer (1,2) and not on 
(3,4). Which set of the two-fold degenerate triplet levels 
is optically active depends on the relative orientation of 
Bq and d2 on the dimer where ^ 0. In Fig.jSjBo || d2 
and in Fig.ElBo J- d2. An additional sphtting of ro„i(±) 
and Top(±) by 0.6 cm~-^ seen in Fig. is because Bq is 
out of {ah) plane by 9°. 

C. Static and dynamic DM in SrCu2 (303)2 

We have shown that the first triplet state energy spec- 
tra are well described with two static DM interactions, di 
and d2. The information about di and d2 is contained 
in the position of energy levels and in the FIR absorption 
line intensities. The inter-dimer di determines the mag- 
netic field dependence of intensities and the triplet state 
level energy splitting. The intra-dimer d2 determines the 
extent of the avoided crossing with iS'2 and the magnetic 
field dependence of intensities near the avoided crossing 
points. Over the magnetic field range of our experiment 
the intensities of the singlet-triplet absorption lines do 
not depend on the dynamic part of the DM interaction, 
because the phonon energies are large compared to the 
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triplet state energy. 

Other inter- and intra-dimer DM interaction compo- 
nents besides di and d2 have been considered to describe 
experimental data>i2iSS These are the in-plane compo- 
nent of the inter-dimer DM d^y and the symmetry- 
forbidden intra-dimer DM in the c-direction. We in- 
cluded dxy and in the two dimer model and found 
that calculations with non-zero dxy and dz give results 
contradicting with the experiment. Our argument, which 
is independent of whether a particular infrared transition 
is allowed or forbidden, relies on the observed and calcu- 
lated crossing - anti-crossing effects between the triplet 
and the singlet states. 

If Bq II c and d^y ^ then rom(+) would have an 
avoided crossing with 5*2 contradicting the experiment, 
where Top{+) anti-crosses with the singlet [Fig.E^a)]. 
Also dz does not give any anti-crossing between S2 and 
Tom{+) or Top(+). In high field nonzero d2 creates an 
avoided crossing between the ground state 6*0 and the 
triplet branch Tom(— ) as observed in the experiment^'* 
while nonzero dxy or dz do not create an avoided cross- 
ing between Sq and T'om(— ) or Top(— ). However, the 
two dimer model does not predict the experimentally 
observedi^ avoided crossing between 5'o and Top{—). 

In Bq II a field orientation both d^y and add, in ad- 
dition to d2, to the avoided crossing of one of the triplet 
components with S2 ■ The experimental data can be fitted 
with a single value ^2 = 1.8 cm^* in both field orienta- 
tions. Bo II a and Bq || c. If da;^ and d^ were comparable 
in magnitude to d2 , then the extent of avoided crossing 
would be different in Bq || a and Bq || c field orientations. 

Our conclusion is that the dominant DM interactions 
are di = 0.7 cm~* and d2 = 1.8 cm~^. In the magnetiza- 
tion plateau state the lattice parameters of SrCu2 (603)2 
may change due to spin-phonon coupling. Our calcula- 
tion of energy levels did not account for static lattice dis- 
tortions and therefore we cannot make any conclusions 
about dxy and d^ and the strength of di and d2 in high 
magnetic fields. 

The intensity of the FIR singlet-triplet transitions de- 
pend on the strength of the dynamic DM and on the fre- 
quency and the oscillator strength of the phonon. Since 
the particular phonons involved in the dynamic DM ef- 
fect in SrCu2 (603)2 are not known we can give only the 
relative strength of dynamic DM interactions. The a- and 
c-polarized singlet-triplet transitions have similar oscilla- 
tor strengths. These are 2.0 cm~^ (Ei || a ) and 1.7 cm~-^ 
(El II c ) if we compare the two lower spectra in Fig.[21 
which have been measured on the same sample by chang- 
ing the direction of the light polarization. The ratio of 
the dynamic DM interactions for the two mechanisms 
is d^a/dzc — \/2 X 2.0/1.7 = 1.5 if we assume that a- 
and c-axis phonons have equal frequencies and oscillator 
strengths. The factor 2 accounts for the a-axis phonon 
creating a dynamic DM only on the dimer with its axis 
perpendicular to Ei. 



D. Staggered g-tensor 

The importance of the staggered g-tensor in 
SrCu2 (603)2 was pointed out by Miyahara et ali^. The 
staggered g-tensor exists in SrCu2(B03)2 because of the 
buckling of Cu-O-B planes below 395 K. It mixes singlet 
and triplet states similar to the static DM interaction d2 . 
The strength of the staggered (7-tensor interaction can be 
estimated and we show that its effect on the energy of 
the spin levels is small compared with the effect of d2. 
The Zeeman term Hzs couples singlet and triplet states 
on a single dimer and is proportional to gsfJ-sBa, where 
ds = {g- — g-) sin ip cos (j) (Ref.— ) . The angle (/> « 6° is 
the buckling angle of the Cu-O-B plane^''. The compo- 
nents ffx and ffz of the Cu ion g-tensor are not known but 
we take gx ^ ga = 1-998 and g^ ~ gc = 2.219 (TableD) 
and get gs = 0.023. The staggered term Hzs increases 
linearly with magnetic field. The largest field where the 
anti-crossing between Tq and 5*2 takes place is 5 T. In this 
field the magnitude of the staggered g-tensor term in the 
hamiltonian is 0.05 cm~^, which is much smaller than the 
static intra-dimer DM term d2 — 1.8 cm^*. We conclude 
that the dominant coupling between the singlet and the 
triplet is due to the static DM interaction d2. 



E. States of bound triplets 

Several states besides the one-triplet excitation are 
infrared- active (TableU)). We showed that the two dimer 
model explains well the energies of the one-triplet states 
and transitions to them. In the two dimer model with 
ji = 24cm~^ and 2j2 = 9.8 cm~* we get several two 
triplet states: a singlet, a triplet, and a quintet of two 
bound triplets at 28.4, 38.2, and 57.8 cm~^, respectively. 

SrCu2 (603)2 has two low energy singlet states 6*1 and 
S2 which both anti-cross with triplet state levels (Fig.^ 
and |SJ). In the two dimer model only one singlet of 
bound triplets is possible and the anti-crossing occurs 
only with Top(±) states. In the experiment an anti- 
crossing is observed between S'2 and Top(-I-), Fig.[Sfa). 
The observed anti-crossing between 5*1 and T'om(— ) can- 
not be explained by the two dimer model. In Section FlI CI 
we show that other DM interactions besides d2 are weak 
or absent in SrCu2(603)2 in the studied Bq range, al- 
though they may have proper symmetry to couple Si 
and To„(-). 

The energy of the 38.2 cm^* triplet in the two dimer 
model is in the range where triplets are present in 
SrCu2(603)2. There is a triplet at 37.5 cm~* labelled 
as Ti (TableP). FIR transitions to this state are active 
in El II a polarization, Fig.|21 In Ei || c the transitions 
are weaker (TablePl. The Ti(0) level is FIR active when 
Bq II c and T'i(±) are active when Bq -L c. All this, 
polarization and magnetic field dependence, is consistent 
with the dynamic DM mechanism of the FIR absorption 
where the dynamic DM is along the c axis. The a-axis 
phonon creates a dynamic DM in the direction parallel to 
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the c axis. The intra-dimer dynamic DM interaction d^a 
does not give any transitions to bound states of triplets. 
We considered a possibility that the a-axis phonon modu- 
lates the static inter-dimer di . We found that the pattern 
of dynamic inter-dimer DM vectors with the same sym- 
metry as di [Fig.^b)] gives selection rules that apply to 
the 37.5 cm"-'^ Ti triplet. Transitions to other states are 
forbidden in the first order of this dynamic DM interac- 
tion. The lattice deformation that creates such a pattern 
of dynamic DM vectors is of symmetry and is not an 
optical phonon; in the Ai symmetry mode Cu atoms on 
ji bond move along the bond in antiphase . We conclude 
that the two dimer model is not sufficient to account for 
transitions to states of bound triplets, except to 5*2. 

Quintet states were observed by high field ESR.^^ 
Their extrapolated zero field energies are in the range 
46 - 58cm~^. There are two Ei || a singlet excitations at 
52.3 and 53.5 cm~^ in this range (TableD). The quintet 
{S — 2) has a TO5 = spin level which has the same 
magnetic field dependence of energy as the S = state. 
However, the observed singlets at 52.3 and 53.5 cm~^ are 
not the ms = components of the quintet. If in one Bq 
field orientation the ms = level is infrared-active then 
in the 90° rotated field orientation other levels, ms = ±1 
or ms — ±2, become active. We studied all possible Bq, 
El orientations relative to crystal axes and did not find 
the splitting of the 52.3 and 53.5 cm"^ excitations in the 
magnetic field although they are one to two orders of 
magnitude more intensive than other magnetic excita- 
tions in FIR spectra. 

We assigned the a axis polarized (Ei || a) 43.0, 52.3, 
and 53.5 cm~^ singlet excitations to magnetic excitations 
because of the magnetic field and temperature depen- 
dence of their energy and intensity. Whether they could 
be phonons activated by magnetic interactions needs a 
further study. 

III. CONCLUSIONS 

In SrCu2(B03)2 the ground state is not exactly a prod- 
uct of singlets on dimcrs as in the Shastry-Sutherland 



model, because the intra-dimer Dzyaloshinskii-Moriya 
interaction d2 mixes the ground singlet state with the 
triplet. From the observed anti-crossing between T'op(+) 
and 52 we get ^2 — 1.8 cm"-'^. This is comparable to 
the inter-dimer DM, di = 0.7 cm"-'^, which determines 
the triplet state energy level zero field splitting. Both di 
and d2 determine the magnetic field dependence of the 
absorption line intensities. 

Although magnetic dipole singlet-triplet transitions 
are allowed by d2 , the experimentally observed polariza- 
tion and magnetic field dependencies of absorption line 
intensities are not described by this interaction. Instead, 
singlet-triplet transitions are allowed by the dynamic DM 
mechanism where the electric field component of FIR 
light couples to a non-symmetric phonon, which creates 
the DM interaction. There are two dynamic DM mech- 
anisms in SrCu2 (603)2- In one case the FIR light cou- 
ples to an a-axis phonon and in the other case to a c-axis 
phonon. This is consistent with the calculations of Cepas 
and ZimanS^ who used a two dimer model in the j2 — 
limit. 

The experiment also yielded information about higher 
triplet and singlet excitations. Several of these absorp- 
tion lines are identified as electric dipole transitions. The 
two dimer cluster is too small to describe these transi- 
tions. Also, we had to use renormalized values of ji and 
j2 to calculate the energy levels because the actual spin 
excitations are delocalized over a larger cluster. Obvi- 
ously a bigger cluster is needed for proper calculation 
of magnetic excitations in SrCu2 (603)2. Nevertheless, 
the two dimer model gives us a good description of the 
one-triplet excitation. 



IV. ACKNOWLEDGMENTS 

We thank G.Blumberg and O. Cepas for fruitful dis- 
cussions. This work was supported by the Estonian Sci- 
ence Foundation Grants Nos. 4926, 4927, and 5553. 



* Electronic address: lroomtom@kbfi.eel 

^ H. Kageyama, K. Yoshimura, R. Stern, N. V. Mushnikov, 
K. Onizuka, M. Kato, K. Kosuge, C. P. Slichter, T. Goto, 
and Y. Ueda, Phys. Rev. Lett. 82, 3168 (1999). 

^ K. Kodama, M. Takigawa, M. Horvatic, C. Berthier, 
H. Kageyama, Y. Ueda, S. Miyahara, F. Becca, and 
F. Mila, Science 298, 395 (2002). 

^ H. Kageyama, M. Nishi, N. Aso, K. Onizuka, T. Yosihama, 
K. Nukui, K. Kodama, K. Kakurai, and Y. Ueda, Phys. 
Rev. Lett 84, 5876 (2000). 

* B. S. Shastry and B. Sutherland, Physica 108B, 1069 
(1981). 

^ S. Miyahara and K. Ueda, J. Phys.: Condens. Matter 15, 



R327 (2003). 

C. Knetter, A. Biihier, E. Miiller-Hartmann, and G. S. 
Uhrig, Phys. Rev. Lett. 85, 3958 (2000). 
^ A. Koga and N. Kawakami, Phys. Rev. Lett. 84, 4461 
(2000). 

* W. Zheng, J. Oitmaa, and C. J. Hamer, Phys. Rev. B 65, 
014408 (2001). 

^ C. H. Chung, J. B. Marston, and S. Sachdev, Phys. Rev. 

B 64, 134407 (2001). 
^° A. Lauchli, S. Wessel, and M. Sigrist, Phys. Rev. B 66, 
014401 (2002). 

" T. Munehisa and Y. Munehisa, J. Phys. Soc. Jpn. 73, 340 
(2004). 



7 



O. Cepas, K. Kakurai, L. P. Regnault, T. Ziman, J. P. 
Boucher, N. Aso, M. Nishi, H. Kageyama, and Y. Ueda, 
Phys. Rev. Lett. 87, 167205 (2001). 

H. Nojiri, H. Kageyama, K. Onizuka, Y. Ueda, and M. Mo- 

tokawa, J. Phys. Soc. Jpn. 68, 2906 (1999). 

H. Nojiri, H. Kageyama, K. Onizuka, Y. Ueda, and M. Mo- 

tokawa, J. Phys. Soc. Jpn. 72, 3243 (2003). 

P. Lemmens, M. Grove, M. Fischer, G. Giintherodt, V. N. 

Kotov, H. Kageyama, K. Onizuka, and Y. Ueda, Phys. 

Rev. Lett. 85, 2605 (2000). 

T. Room, U. Nagel, E. Lippmaa, H. Kageyama, 
K. Onizuka, and Y. Ueda, Phys. Rev. B 61, 14342 (2000). 
C. Knetter and G. S. Uhrig, Phys. Rev. Lett. 92, 027204 
(2004). 

A. Zorko, D. Arcon, H. van Tol, L. C. Brunei, and 

H. Kageyama, Phys. Rev. B 69, 174420 (2004). 

S. Miyashita and A. Ogasahara, J. Phys. Soc. J. 72, 2350 

(2003). 

J. Redhammer, 
Lemmens, A. 
F. Hiining, et al 



P. Roussel, G. Heger, 
lonescu, M. Grove, 
Eur. Phys. J. B19, 507 



K. Sparta, G. 
G. Roth, P. 

G. Giintherodt, 
(2001). 

^1 K.-Y. Choi, Y. G. Pashkevich, K. V. Lamonova, 

H. Kageyama, Y. Ueda, and P. Lemmens, Phys. Rev. B 
68, 104418 (2003). 

B. Wolf, S. Zherlitsyn, S. Schmidt, B. Liithi, H. Kageyama, 
and Y. Ueda, Phys. Rev. Lett 86, 4847 (2001). 
O. Cepas and T. Ziman, Phys. Rev. B 70, 024404 (2004). 
H. Kageyama, K. Onizuka, T. Yamauchi, and Y. Ueda, J. 
Crystal Growth 206, 65 (1999). 

T. Room, D. Hiivonen, U. Nagel, Y.-J. Wang, and R. K. 
Kremer, Phys. Rev. B 69, 144410 (2004). 
^® Weak 52.3 and 53.5 cm~^ lines are present in Ei || c spec- 
tra because the crystal was misaligned relative to the light 
k vector; that created a small in-plane Ei component. 
As estimated from the residual intensity of the 52.3 and 
53.5 cm~^ lines in Ei || c polarization the k vector was out 
of the (ab) plane by 12 ± 2° . The same applies to the Bo 
field in Faraday configuration, Bo || k. Misalignement of 
9 ± 1° is calculated from the splitting of the triplet com- 
ponents, Fig.El 

L. Shekhtman, O. Entin-Wohlman, and A. Aharony, Phys. 
Rev. Lett. 69, 836 (1992). 

L. Shekhtman, A. Aharony, and O. Entin-Wohlman, Phys. 

Rev. B 47, 174 (1993). 
^® G. Jorge, R. Stern, M. Jaime, N. Harrison, J. Bonca, S. E. 

Shawish, C. Batista, H. Dabkowska, and B. Gaulin, cond- 

mat/0309534 (2003). 
^° S. Miyahara, F. Mila, K. Kodama, M. Takigawa, M. Hor- 

vatic, C. Berthier, H. Kageyama, and Y. Ueda, J. Phys.: 

Condens. Matter 16, S911 (2004). 




(b) 



Cu2+, S=l/2 t 

. yllb 




FIG. 1: Cluster with two dimers (1,2) and (3,4). (a) Dimer 
(3, 4) and four nearest-neighbor dimers. The thin dashed 
line shows the two dimer cluster boundary. Thin solid lines 
show the distortion of Cu-Cu superexchange bonds due to the 
buckling of Cu-O-B planes. Thick solid and dashed lines are 
the inter- and intra-dimer superexchange constants ji and j2; 
inter-dimer DM vectors (di, solid arrow) are in the c direc- 
tion and intra-dimer DM vectors (d2, empty arrow) in the 
(ab) plane along a and 6 axis, (b) The two dimer model after 
the periodic boundary condition has been applied; inter-dimer 
interactions have doubled. 
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FIG. 2; Differential absorption in Ei ± c (two upper curves) 
and El ||c (lower curve) polarization. Spectra have been off- 
set in vertical direction. 
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FIG. 3: Differential absorption spectra in magnetic field Bo || 
a at 4.4 K. Vertical offset equals to the magnetic field value 
in Tesla. 




FIG. 4: (color online). Magnetic field dependence of line positions and line areas in Ei ||c polarization at 4.4K; (a), (b) Bo || 
(c), (d) Bo II a. Solid lines are the results of the calculation based on the two dimer model: ji =24cm~^, 2j2 =9.8cm^ 
2di =1.4 cm" ^, and d2 =1.8 cm" ^. Dashed lines in panels (a) and (c) are fits with parameters given in TableQ] 
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FIG. 5: (color online). Magnetic field dependence of line positions and line areas in Ei ||a polarization at 4.4K; (a), (b) Bo ||c; 
(c), (d) Bo II a. Solid lines are the results of the calculation based on the two dimer model: ji =24cm~^, 2j2 =9.8cm~^, 
2di =1.4 cm~ , and d2 =1.8 cm~ . Dashed lines in panels (a) and (c) are fits with parameters given in TableEl The sohd line 
in panel (b) is the sum of two theoretical line areas of So to rom(O) and to rop(O) transitions shown by dashed lines. Dashed 
lines in (d) are eye guides (see text). The 18 T point in panels (a) and (b) was measured at 1.8 K. 
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FIG. 6: (color online). Line positions (a) and line areas (b) 
in El II a and Bo || b configuration at 4.4K. The lines are 
results of the calculation based on the two dimer model and 
dynamic DM interaction. The additional splitting of triplet 
components (triangles) is caused by the magnetic field Bo 
being misaligned by 9° out of the (ab) plane. In panel (b) 
the line area (triangles up or triangles down) is a sum of line 
areas of split components. 




FIG. 7: Intra-dimer dynamic DM interactions. A lattice dis- 
tortion with the normal coordinate Q (solid arrow) creates 
an intra-dimer DM interaction da (empty arrow). The c-axis 
phonon creates a dynamic DM interaction on both dimers 
while the o-axis phonon affects the dimer (1,2) only. 
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